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We study p-wave pairing in a two-component Fermi system with unequal population across weak- 
coupling BCS to strong-coupling BEC regimes. We find a rich rris = spin triplet p-wave superfluid 
ground state structure as a function of population imbalance. Under a phase stability condition, the 
"global" energy minimum is given by a multitude of "mixed" SF states formed of linear combinations 
of m = ±1,0 sub-states of the t= \ orbital angular momentum state. Except for the "pure" SF 
states, (i = \,m = ±1), other states exhibit oscillation in energy with the relative phase between 
the constituent gap amplitudes. We also find states with "local" energy minimum that can be stable 
at higher polarizations, suggesting a quantum phase transition between the "global" and "local" 
minima phases driven by polarization. The local and global minimum states may be associated with 
Morse and non-Morse critical points. We discuss possible consequences for experiments. 
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Introduction: Over past several years there has been 
sustained interest in paired fermion ground states with 
unconventional pairing symmetry. Among these are p- 
wave spin triplet condensates. Correlated electron sys- 
tems, such as, SrRu204 [1] and ferromagnetic supercon- 
ductors P] are believed to possess p-wave triplet symme- 
try. Fermi systems with unequal species population, as 
quark matter [3j, magnetic field induced organic super- 
conductors [3] and cold fermion systems with unequally 
populated hyperfine states [5], add a fascinating dimen- 
sion. Discovery of s-wave superfluidity in cold atoms [6i 
subjected to s-wave Feshbach resonance (FR), and obser- 
vations of p-wave FR [THIO] in ^Li and ^^K had raised 
the prospect for observing p-wave superfluidity in cold 
Fermi gases. While this may be difficult to achieve, al- 
ternate methods [TT] and optical lattices [T^ [13] offer 
encouraging prospects. There has been theory work on 
p-wave superfluidity in BEC-BCS crossover region for a 
single-component Fermi gas |14H17j , and two-component 
case jl8] with equal-species population. However, the- 
oretical work on p-wave pairing for unequal population 
has been limited [131 |T3] . 

In this paper, we conduct analytical and numerical 
study of p-wave superfluidity in a two-component Fermi 
system with unequal population, across weak-coupling 
BCS to strong coupling BEC regimes. We consider the 
case where inter-species pairing interaction is dominant. 
While this is interesting to study on general grounds, 
we are also motivated by work |13j on population- 
imbalanced Fermi mixtures in optical lattices, and the 
observation [HI dS] that unlike liquid '^He, pairing in- 
teraction in cold atoms may be highly anisotropic in 
"spin" -space, ("spin" referring to hyperfine states). For 
example, in ^Li, when the hyperfine pair |ms,ml, >= 
|l/2, — 1/2 > is at resonance, the pairs |l/2,l/2 > and 
I — 1/2,-1/2 > would not be. So, intra-spin interac- 
tions can be neglected. Pairs in p-wave superfluids with 



unlike "spin" components can, however, have different 
£ = 1 components: m = ±1,0. Then the gap parameters 
'-^e=i,m are related to the spherical harmonics, Yi,rn=±i,o- 

Our key results are: We find a rich superfluid (SF) 
ground state (GS) structure. Under a condition on 
relative phase between pairing amplitudes A„i=o,±i 
{Ag^i^rn = A,„), a multitude of "mixed" superfluid states 
of the form uqAq ± aiAi ± a_iA_i are found to be de- 
generate with A-ti; these give the global minimum GS 
energy. Additionally, we find states with local minimum 
energy. We provide a geometric representation of the 
states (Fig. 1). Our T=0 polarization (P) vs p-wave 
coupling phase diagram (Fig. 2) shows two SF phases, 
comprising respectively of states with global and local 
minimum energy, and a region of phase separation. This 
suggests the intriguing possibility of a quantum phase 
transition between the two SF phases, driven by polar- 
ization. The P GS structure is preserved in the 
P — > limit; hence richer than previous P = results 
[18) . The energies of the "mixed" states oscillate with 
the relative phase angle (Fig. 3); this may be interesting 
to explore experimentally. 

Model: We consider a system of fermions with unequal 
"spin"((t,4.) = (1/2,-1/2)) population, but with equal 
mass. Pairing interactions in all ^ = 1 channels (to — 
0, ±1) are assumed equal; intra-species interactions (^^ 
or J,J,) are set to zero. Since we adjust self-consistently 
the chemical potential with the strength and sign of the 
coupling, in our fermion-only model, molecules appear 
naturally as 2-fermion bound states. The 5 = 1, — 
triplet pairing Hamiltonian is given by: 

kcr 

+ ^kk'4+q/2t^-k+q/2i'^-k'+q/2;Ck'+q/2t (1) 

kk'q 

where Cko- (cj^g.) the annihilation (creation) opera- 
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tor for a fermion with momentum k, kinetic energy 
Ckcr=ek;— Mcr, Spin u] /icr is the species chemical potential, 
and ek=fi-^fc^/2m. Vkk' is the pairing interaction. 

We consider condensate pairs with zero center-of-mass 
momentum (q=0). While non-zero q is interesting from 
the perspective of FFLO [20] states, the problem of q = 
unconventional pairing in unequal-population systems is 
rich in itself. T-L is mean-field (MF) decoupled via S — 
l,ms = "spin" -triplet pairing gap function A^^(k) = 
^(k)=-Ek' ^kk' (c-k'iCk't) giving: 

H^^ = E^k.cLck. - 5: [A(k)4,cl,,+H.c." 

ko- k 

- ^kk' < 4^clk^ X c_k'iCk't > (2) 

k,k' 

To obtain GS energy of the MF Hamiltonian ^ we 
use the equation-of-motion method jSl] |22] for the 
imaginary-time (r = it) normal {Gaa') and anomalous 
(i^o-cr') Green's functions [53]: 

+ A_,,(k)^^_,,,(k,T), (3) 
drF„a'{KT) = C-k<T-F<^CT'(k,r) 

+ A;_,(k)G_..,(k,r). (4) 

Eqs. ^ and Q are Fourier transformed with 
T^iujn=v and dr^—v, where iujn—{2n-\'l)-KT are 
Matsubara frequencies. This gives G^crlk, i^) — 
-(^^k-a+i^)/ [(ek^-i/)(e-k-^+i') + |A^,_<,(k)|2] , and 
F,_,(k,z.)=A*_,(k)/[(ek-.-;^)(e-k.+^^) + |A,_,(k)|2] 
= -F_,,(k, v). 

In terms of radial and angular parts, Vkk' = 
(W3)E„,^^,fc'n.™(k)yi:™(k') (k=(0»). Using 
(c-k'^Ck'^) = -F^.j.(k', r=0 ) and Fourier transforming, 

we obtain, A(k)=-(l/fcBT) E^^, Vkk'^^^lk', i^)e''0+ . 
We take Vj.^' — Xw{k)w{k'), a separable form 
chosen for convenience [121 dHl [211; ^his docs not 
change physics qualitatively. Then the gap function: 
A(k) = J2m^mw{k)YiM = wik)A{k), where A,„ = 

-(l/fcsT) E.k' A«^(fc')i"i:„(k')nT(^'' Obtaining 
F^^(fc, i/) and Gcro- from Eqs. (3)-(4), we get for the 
£=1;to = 0,±1 gap amplitudes, 

A™ = -Xj2y^ik)YUmmnF{E^ynF{E+)]/{2E,) 

k 

(5) 

and for the particle number densities, 

A; fc 

Above, h = (Cfe; - efet)/2 = (Mt - M;)/2, - 
K2 + |A2(k)|]i/2, with efe = (efeT+CH)/2=efe - /i and 
/i=(/x^-|-/i^)/2; np{E^) are Fermi functions with = 



—h ± At T=0 and for /i>0, the expectation value of 
the grand canonical potential is given by: 

+ J2 m-h)9{-E,+h)} - 1 ^ \A„f (7) 

k ^ rn 

Following standard practice, coupling A has been ex- 
pressed in terms of the "regularized" interaction g: 
l/g=l/A-h(l/(27rfi)3) / w^{k)d^k/2ek [25|. For w{k), we 
adopt the N-SR form w{k)=kok/{kl + P). 

Ground State Analysis: We assume < |A(k)| <C ^k, 
and expand E^ in Eq.Q in powers of |A(k)| to 4th 
order: E^ ^ |a|[l+|A(k)|V2|ap-|A(k)|V8|a|^]. 
Since £,k ^ |A(k)|, Cfe = i^k = m) is excluded 
from the expansion [^B]. Substituting this in Eq.Q 
and converting momentum sums to integrals we 
obtain: Eg = aJ\A{k)\^dn + /3 J \A{k)\'^dn + 7, 
where a ^ {l/{2Trh)^) Jj^^i_^w'^{k)k'^dk/2ek + 
{l/{27Thr) wHk)k'dk{l/2ek-l/2\^k\) ~ 1/5, 

/3 = {l/{2TThr)J^^^^wHk)k^dk/8\^k\' and 
7 = (47r/(27r/i)^) Jp^^,^i£.k—h)k'^dk. Existence of a 
superfluid phase requires a<0, otherwise minimization 
of Eg forces the gap to vanish. The polarization- 
dependence {P = {n^ — ni) / {n^ + n^)) of Eg is contained 
in a,/3,^, which depend on h, fi and w{k); a alone 
depends explicitly on coupling g. Accordingly, the A,„'s 
depend on P. 

For fixed (/i, h), and Ek to 0(A(k)''), we obtain an 
analytic expression for the ground state energy, Eq, in 
terms of the pairing amplitudes Am- 

Eg ^-^+1 + 213 (^|AoP + \A,f + |A_ip + 

+/3(|Aop-2|Ai||A_i|)2+4/3(l-t)|AonAi||A_i|(8) 

where t=cos6', with 9 = <(AoAJ, A;5A_i) = 
2(/)o— 01— (/)m's are phases associated with the gap 
amplitudes A^. Minimizing Eq in ^ with respect to 
Am's, we find that for a stable p-wave superfluid phase 
to exist, the following conditions have to be satisfied si- 
multaneously: 

(a) |Ao|2-f|Ai|2+|A_i|2 = -a/4/3 = i?2 (sphere) 

(b) |Ao|2-2|Ai||A_i| =0 (plane) (9) 

(c) (l-t)|Aop|Ai||A_i| ^0 

For all three A^'s non-zero, conditions (a) and (b) 
give a semicircle formed by the intersection of the surface 
of a sphere of radius R, [R^ — —a/4/3) in Ai, A_i, Aq 
space, with a plane defined by |Ai| -I- |A_i|=i?. Points 
spanning the semicircle represent a multitude of " mixed" 
superfluid states of the form CoAq + aiAi + a_iA_i; 
ao,ai,a_i being constants. A 3D geometric represen- 
tation of the states is shown in Fig.l. Condition (c) 
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FIG. 1; (Color online) Geometric representation of p-wave 
states. Semicircle formed by intersection of sphere surface 
with plane represent a continuum of states exhibiting global 
energy minimum (see text). Plane containing semicircle is 
also shown (green). O (on Ao axis) and M (in Ai-A_i plane) 
are states with local energy minimum. Am(»7i ~ 0, ±l)'s nor- 
malized to sphere radius \J —aj^^fi; a < 0. 

imposes the constraint, t = cos6 = 1, i.e. the relative 
phase 9 ~ (</>o— (/ii) — (<?!)-i— 0o) = 2mr. t=l corresponds 
to a parallel orientation of vectors AqA^ and AqA_i. 
The ground state (GS) global minimum energy, 

E^(j = ~a^/8l3 + 7, (a < 0), (10) 

is completely determined by a, /3, 7 given earlier. 

The phase condition (c) is always satisfied for the 
states, A±i, at the endpoints A,B of the semicircle. But, 
energies of other states on the semicircle oscillate with 
6, and attains only for t = 1, i.e. 9 = 2mT (Fig. 
3a). The oscillation amplitudes depends on the specific 
state; state C at the top of the projected figure (Fig.l) 
has the maximum amplitude. Representative states 
A,B,C on semicircle are given by: (A) |Ao| = |A_i|=0 and 
|Ai|2=-a/4/3; (B) |Ao| = |AihO and | A_i |2=-a/4/3; 
(C) |Aop=-a/8/3 and |Aip=|A_i|2= - a/16;9. 

We obtain other minimum energy solutions by consid- 
ering special cases in Eq. [sj (i) |Ai| = = |A_i|, 
|Ao| 7^ 0: we find a unique solution, |Aop — —a/6/3; 
point O in Fig.l. (ii) |Ao| = 0: this gives equations for 
two ellipses in Ai — A_i plane: 

|Ai|2 + 2|A_i|2 + ^ =0 (ellipsel) 

2|Ai|2 + |A_i|2 + ^=0 (ellipse2) (11) 

Intersection of the two ellipses gives the solution: | Ai p — 
|A_i|2 = -a/12/3, |Ao| = 0; point M in Fig.l. Like 
A and B, states O and M automatically satisfy phase 
constraint (Eq. 9(c)). These exhaust all possible inde- 
pendent solutions. Incidentally, states A (B) respectively 



lie at the intersection of ellipse 1 (ellipse 2), and sphere 
and plane of Fig.l. The energies for O and M are equal, 
and lie higher than Eq ; we denote this as local minimum 
energy: 

E'"" = -a'^/Up + 7 (a< 0) (12) 

We checked the stability of the states by examining 
the stability matrix or Hessian, H: {d^EQ/dArmd^m )- 
det[H] = for states A, B, and det[H] > for states O, 
M; in the language of Catastrophe Theory [27], A, B can 
be identified with non-Morse "critical" points, and O, M 
with Morse "critical" points ,28, . 

As can be seen from the expression for 7, and the defi- 
nition of /i , /i, in the limit of zero polarization (P — 0), 
the GS energy coefhcient 7 — > 0, while a < 0, /3 > 
attain different numerical values. Accordingly, expres- 
sions for E'§ and E^"" remain essentially the same for 
P — > 0, but with different numerical values. States A 
and B (equivalent to A by symmetry) correspond to the 
finding of Ref. [H] for P = 0. 

Numerical Calculations: Detailed numerical study of 
the p-wave superfluid states confirms our analytical re- 
sults, and reveals additional features. We solve self- 
consistently three gap equations (Eq. (5)) and two num- 
ber equations (Eq. (6)), for fixed population imbalance 
P = (n-t- — n\^)/{n^ + n\^), and p-wave coupling param- 
eter g. (g is related to the triplet scattering parameter: 
g = 25kpat/8TT for a cutoff ko = lOkp)- This gives us 
the gap amplitudes A„i, and the chemical potentials ^a- 
for different P and l/kpUt- Using these, we obtain the 
GS energy Eg from Eq. (7), and the coefficients a,/3,7 
that determine Eg] agreement between our numerical and 
analytical results for Eg is good. We note that fixed 
(P, n — + n^) is equivalent to fixed (/i, h) used in an- 
alytic study above. 

To check for stability of the p-wave states obtained 
from numerical solutions, we enforce that the stabil- 
ity matrix {d^ Eq / dA^idAjn-) is positive definite; and 
6p/Sh > 0. Based on this, we construct [55] a polar- 
ization (P) - coupling (l/kpat) phase diagram, Fig. 2, 
in BEC-BCS regimes. SFl denotes a stable superfluid 
phase, corresponding to the states on the 'semicircle' 
(Fig.l) that give GS global minimum. Eg. With in- 
creased P, polarization-dependent parameters, a,/3,7, 
that determine Eq and Eq'^, change so that the super- 
fluid phase SF2, corresponding to states with local min- 
imum, attain energy lower than that of SFl, and hence 
becomes stable. This suggests the interesting possibility 
that at T=0, polarization may drive a quantum phase 
transition from SFl to SF2. It may also be possible to 
access SF2 at T^O. At even larger polarizations, SF2 be- 
comes unstable to phase separation (PS). SFl, SF2, and 
PS occupy a relatively much narrower part of the phase 
diagram on the BCS side compared to the BEG side. In 
our two-component system with inter-species interaction, 
PS persists into full polarization, P=l. This is because at 
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FIG. 2: (Color online) Polarization P vs p-wave coupling 
— l/{kpat) phase diagram, showing normal (N), p-wave su- 
perfluid phases (SFl; SF2), phase separation (PS). Unitarity 
limit is shown by the dashed vertical line. Inset: Calculated 
chemical potential vs coupling across BCS/BEC regimes for 
P = 0.1 (solid line), 0.5 (dashed line). Fermi energy, ep, is 
given by t-i/ep ~ 1 line. 



P=l, the system is essentially a one-component system; 
absence of minority-species fermions makes inter-species 
interaction inoperative. Such a system can exhibit su- 
perfluidity only with intra-species interactions. 

The inset in Fig. 2 shows the calculated behavior of 
chemical potential fi across BEC-BCS regimes. It devi- 
ates significantly from the Fermi energy {fJ-fep — 1 hne) 
in a wider region around the BEC-BCS crossover, even 
on the BCS side, and drops much more rapidly to nega- 
tive values on the BEC side compared to the s-wave case. 
For sufficiently weak coupling (BCS regime), fi ep. 

Except for A, B, the continuum of "mixed" SF states 
on the semicircle (Fig.l) are characterized by relative 
phase 2mT < < {2n + 2)tt, so that their energies oscil- 
late with 6. Numerical results for three cases are shown 
in Fig 3(a). For 9 ^ 2mr, energies lie higher than the 
GS global minimum; maximum amplitude occurring for 
state C in Fig.l. This raises the possibility of observing 
^-oscillation in the states on the semicircle using phase 
sensitive experimental technique(s). They could also be 
accessed at finite-T. Figs. 3b, 3c respectively show the 
variation of Eg with coupling g at fixed polarizations P, 
and with P for fixed g. Eg is normalized to the Fermi en- 
ergy ep = h'^kp/2m, with 2kp = kp^ + kp^. For a given 
P, Eg becomes less negative as g approaches unitarity; 
the trend is more noticeable for smaller P's. The energy 
curves in Fig. 3(b) would terminate at some non-zero 
coupling, and hence not expected to cross. Extrapolated 
to normal regime, energies corresponding to different po- 
larizations would be different, and also not cross. For a 
given g, Eg lies higher for smaller P, presumably due to 
the lower majority-species band becoming progressively 
more occupied with increasing P, thereby lowering Eg 



FIG. 3: (Color online) (a) Calculated ground state energy 
(scaled to global minimum) vs. relative phase angle 9 
(see text). Curve with maximal amplitude corresponds to 
state C in Fig 1: |Aop=-a/8^, j Aip=j A_ip=-a/16^. 
Other two curves correspond to other representative states 
on the arc in Fig. 1: | Ao|^=-a/11.3/3, | Ai |^=-q/76.3/3, 
A_i|2=-a/6.73/3; | Ao|^=-q/16/3, | Ai|2=-q/186/3, 
I A_i|^=— a/5.49/3, (b) Calculated ground state energy Eg 
vs. coupling g for different polarizations P. (c) Calculated 
Eg vs. P for different g. 



with increasing polarization. For small P, Eg becomes 
less negative with increasing g. This trend is reversed 
for P > 0.3. The crossing at P ss 0.3 is suggestive of a 
possible scaling using stretches in P and Eg. 

Summary: With inter-species pairing interaction in 
a population-imbalanced Fermi system, we find a rich 
p-wave superfluid GS structure involving various sub- 
states of orbital angular momentum i = 1. States giv- 
ing global and local energy minimum are associated with 
non-Morse and Morse critical points. The 3D geometric 
depiction adds insight into the superfiuid states. Our nu- 
merical calculations suggest a quantum phase transition 
between two SF phases driven by polarization. Energies 
of a multitude of "mixed" SF states show oscillations 
with a relative phase angle, that could be observable in 
phase-sensitive experiments. Insight into the nature of 
the orbital part of our superfluid states may be gained 
from measuring angular dependence of momentum dis- 
tributions; molecular spectroscopy using light radiation; 
or possibly zero sound attenuation. Our work should 
apply to other unequal-population Fermi systems, espe- 
cially where p-wave intra-species couplings are small or 
negligible. 
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